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Abstract 

Relying on the collinear factorization approach, we demonstrate that HI and ZEUS measurements 
of exclusive light vector meson and photon electroproduction cross sections can be simultaneously 
described for photon virtualities of Q > 2 GeV. Our findings reveal that quark exchanges are 
important in this small Xbj region and that in leading order approximation the gluonic skewness 
ratio is much smaller than one. 
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1 Introduction 



The HI and ZEUS collaborations intensively studied exclusive electroproduction reactions, such 
as deeply virtual meson production (DVMP) ofp° [Il|2l[3ll2|5l[6l[7], [HIEllIOllIIlE], uj [I2] 
and J/i^ [2l Ha m [15] and deeply virtual Compton scattering (DVCS) [lEl [HI [TBI HH], in the 
small XBj kinematics. Phenomenologically, the DVMP processes have been widely discussed with 
color dipole models, see e.g., Refs. [201 1211 [221 ISSl [211 [251 126], where a rather good description of 
measurements was reported [3 [6]. The underlying idea of such models, applicable at small xbj, 
is that the virtual photon splits into a quark-antiquark pair that interacts via a gluonic t-channel 
exchange with the proton [271 [2H1 [221 [201 [31] • On the other hand a coUinear factorization theorem 
was elaborated that allows to resolve the partonic content in deeply virtual processes by means 
of perturbation theory [32j. It states that for a longitudinally polarized photon exchange the 
DVMP amplitudes factorize into mesonic distribution amplitudes (DAs) and generalized parton 
distributions (GPDs), which are convoluted with a partonic coefficient. Furthermore, the partonic 
part, including the changes of GPDs and DAs under scale variation, can be systematically evalu- 
ated in perturbation theory and is nowadays known at next-to-leading order (NLO) [33] [311 [35] . 
Contrarily to the color dipole model, in the coUinear factorization approach both quark and gluon 
t-channel exchanges are important for the experimentally accessible small xbj region. Also a GPD 
inspired hand-bag model for the DVMP description has been proposed by Goloskokov and Kroll 
[361 [37]. Here, the proton to proton transition is described by a coUinear GPD while the parton to 
meson transition part includes also transverse degrees of freedom. Although much attention has 
been given to HI and ZEUS measurements within the aforementioned models, the description of 
these DVMP data has so far not been explored in the coUinear factorization approach. 

In the coUinear factorization approach the t-dependencies of the longitudinally DVMP and 
DVCS cross sections arise from those of GPDs. Hence, if these processes would be dominated by 
gluon exchanges in the small xbj region, the t- dependencies of the various cross sections should 
become universal. Experimentally, the exponential t-slope of DVMP cross sections for light vector 
mesons decreases with growing and approaches at moderate photon virtuality the DVCS one, 
however, they are with B{Q^ ^ A GeV^) ~ 6/GeV^ larger than the B^/"^ ~ 4.5/GeV2 slope of J/^ 
electroproduction, see, e.g.. Fig. 4 in [38]. Taking the universal t-dependency criteria literally, it 
has been argued from the experimental findings that the onset of the perturbative regime appears 
at rather large photon virtuality of ~ 15 GeV^ or so. Theoretically, this is somehow supported 
by numerical studies in which model dependent NLO corrections turn out to be large [311 IM] and, 
moreover, at this order the residual factorization and renormalization scale dependencies might 
be still rather strong. Note that these scale setting uncertainties should be maximal at leading 
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|. On the other hand the DVCS amphtude is in the colhnear factorization approach 



dominated by quark exchanges rather a gluonic one and the cross section measurements can be 
well described at LO and beyond [ID], where radiative corrections can be considered as moderate 
|¥T] . In turn we might argue that DVMP of light vector mesons in the small xbj region can 
be perturbatively described already for a photon virtuality of > 4GeV^, where the t-slope 
deviations of light and heavy meson vector electroproduction cross sections might be attributed 
to differences in the transverse distribution of sea quarks and gluons. 

The most straightforward method to judge on the perturbative description of these processes 
is a global fit to all of them. For doing so, one needs a flexible GPD parameterization, which is 
elaborated in terms of a Mellin-Barnes integral transformation that maps conformal GPD mo- 
ments into the momentum fraction space [HI HI]. Although the NLO corrections are known in 
this conformal representation B3] the software tools for such DVMP fits are presently under 
development. To get a first insight in the phenomenological description of DVMP processes in the 
small XBj region by means of the collinear factorization approach, we restrict ourselves to the LO 
approximation and fit flexible GPD models to experimental measurements. 

The remainder of the article is as follows. In Sect. [2] we shortly introduce the theoretical 
formalism and set up our GPD models in terms of conformal moments. In Sect. [3] we confront 
then this GPD framework with DVMP of p° and as well as DVCS measurements from the HI 
[SI HZl HSj and ZEUS [IH [3 HSl HB] collaborations. We give predictions from DVCS fits P] and 
the hand-bag model [37] , and confront them in return with DVMP and DVCS measurements. We 
also present the first simultaneous GPD fits to DVCS and DVMP measurements, which illustrate 
that in spite of various theoretical and experimental uncertainties the collinear framework might 
be applicable in the small xbj region for > 4GeV^. Finally, we summarize and conclude. 

2 Formalism and GPD modeling 

The DVCS cross section at small xbj is dominated by the target helicity conserved CFF Ti: 



where a ~ 1/137 is the electromagnetic fine structure constant and the ellipse stays for kine- 
matically suppressed contributions, which include besides t/4M^ and xbj proportional terms also 
non-dominant twist-two CFFs and S, see, e.g., Ref. [ID]. To LO accuracy in the running cou- 
pling constant the CFF 7/ is decomposed in terms of charge even partonic CFFs, which we 

^At this order, e.g., the ambiguous setting of the factorization scale fj, in the DA and GPD is not compensated 
by a change of the hard-scattering coefficient, see below ([5]). 
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denote in the following as "CFF" . For four active quarks we write 
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where the "CFFs" arise from the convolution of the corresponding GPDs with the LO coefficient, 
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Here, we express the scaling variable ^ = a;Bj/(2 — xbj) by the Bjorken variable. 

The longitudinally polarized DVMP cross section reads in the kinematics of interest as follows: 



da^lP^^P Tw-2 . x^ 



dt 



(4) 



where [K^^ is a helicity conserved transition form factor (TFF) and the kinematically suppressed 
contributions, indicated as ellipse, include also the target spin-flip TFF |l5l [Ml |39]. The 
dominant TFF CK^^ factorizes at leading twist-two and at LO accuracy in a^. 
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into the inverse moment of the vector meson DA 
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and the Ti^^ amplitude that contains the GPDs. Note that the TFF ([5]) is proportional to ag and 
to fv/Q- Hence, in this approximation the residual renormalization scale fir and factorization 
scale fi dependencies are of order a^. For light neutral vector mesons these T-L^^ amplitudes are 
decomposed as follows 
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where the charge even quark "CFFs" are given in ([3]) and the gluonic one is deflned as 



'H'^{xB,t,^ 
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The prefactors in (I7H9]) arise from both the electrical charges of quarks and the flavor content 
of the meson DA. Furthermore, we take in ([5]) for the meson decay constants fpO ^ 209 MeV, 
^ ^ 195 MeV, ^ 221 MeV. 
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To conveniently treat evolution, it is rather popular to expand the DA in terms of conformal 
partial waves (CPWs). This provides us for the inverse moment (jS]) the series 
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^""{f^') = J2Ek{fi',fil)vl{fil) with 

even 

in terms of CPW amplitudes <f^. Their scale dependency is governed by the evolution operator 



( ln(^VA 
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which is defined in terms of the LO anomalous dimensions 7^°'', where S\{}z) is the first order 
harmonic sum, the QCD scale parameter, which foinf = 4 active quarks is set to Aqcd = 218 MeV, 
and the renormalization coefficient /3o = 2n//3 — 11 = —25/3 of the running coupling. The inverse 
moment flTT]) is normalized to one in the asymptotic limit /i^ — t- 00. In the following we adopt this 
asymptotic value also at the input scale. Note, however, that sum-rule estimates [IH] provide for 
the p meson DA the second conformal moment 

a^(/i = IGeV) = 0.18 ±0.1, 

indicating a moderate deviation from the asymptotic value X = 1. 

We employ also the CPW expansion for "CFFs" and GPDs 1^ ITT]. However, in contrast 
to the inverse DA moment ([6]), given as a series (fTTj) . "CFFs" fl3l[T0|) are now represented in terms 
of a Mellin-Barnes integral, which reads in the flavor non-singlet sector as 
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r(|)r(3 + j) 

is the hard-scattering amplitude at LO in the charge even sector, the evolution operator Ej is 
defined in ([12]), and Hf^^'^\^,t, /ip) are the conformal GPD moments, analytically continued from 
the odd ones j = 1, 3, ■ ■ ■ . The GPD moments are specified by partial wave amplitudes Hj j that 
appear in their S0(3) t-channel partial wave expansion [18], 

i+i 

if,(r/,t,/i2) = 5^r/^+i-^iJ,-j(t,^2)^V(^) j = l,3,---. 



J-0 
even 



where d'^{r]) are (some) Wigner rotation matrices, labeled by t-channel angular momentum J 
and normalized to one in the limit — > 0. In the forward case, A = 0, the leading SO (3) PW 
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amplitudes, i.e., J = j + 1, are constrained by the Mellin moments of common parton distribution 
functions (PDFs), 

•1 



Hj{r] = 0,t = 0, fi^) = Hj^j+i{t = 0, fi' 



dxx^q{x, fi ) . 



In our GPD model we implement the skewness effect by taking three SO (3) partial waves, where 
the two non-leading ones are expressed by the leading one, multiplied by the strength S2k- 

Hj+2kj+i{t,l^^) =S2kHj{t,lJ^), A; = {0,1, 2}, sq = 1 . (15) 

We emphasize that this effective model allows us to control both the normalization of the "OFF" 
and its change under evolution. Finally, interchanging the J summation with the integration over 
j provides the formula 

2 -I <"c+ioo " / 

dj i + tan ( ■ ^ 



k=0 



Cj+2k Ej^2k{fJ' ,/^o) S2k Hj{t, fiQ 



(16) 



which is used for numerical evaluation. Note that we neglected here the skewness dependency of 
Wigner's rotation matrices, appearing in the CFF ([T51) . which is a save procedure in the small 
region [H]. 

In the flavor singlet sector the quark combination 



and gluon 'H'^ "CFFs" will mix with each other. Adopting the conventions of Ref. 
introduce two dimensional vectors for "CFFs" and GPD moments 



we 



In the forward case the moments 

Hf 
Hf 



and Hi 
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are normalized to the Mellin moments of the common flavor singlet quark (S) and gluon (G) 
PDFs. Analogously to f|T6|) . we model the small ^ behavior of the singlet "CFF" as following 
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where ^ denotes matrix multiplication. Here, the LO coefficient matrix reads 

Lo 2^+ir(j + 5/2) / 1 



r(3/2)r(j + 3) \ 



(19) 



the evolution operator Kj is specified as a two dimensional matrix in Sect. 4.2 of Ref. [H], and 
the two parameters of the diagonal matrix 

= I ^ (20) 

control the skewness effects of the singlet quark and gluon GPDs. 

To have a simple comparison with findings from DVCS fits [H], we choose as input scale 
^ 4 GeV^ and equate renormalization and factorization scales with the photon virtuality. For 
charge even contributions, i.e., nonnegative integer j is odd, the (conformal) moments contain 
contributions from the valence-, sea-, and anti-quarks. Supposing that the sea quark GPD has 
the same functional dependence as the anti-quark ones, we write in analogy to PDF terminology 

^ + 2i/| ^ 2Hl (21) 

where in the small xbj region we can safely neglect valence contributions. Furthermore, we simply 
assume that the functional form is fiavor independent, and hence fl2T|) can be expressed by the 
total sea contribution 

Hf^^ ^ Smp, Hp = 2 J2 S^ = l, (22) 

q=u,d,--- q=u,d,--- 

where S'^ is the fiavor asymmetry parameter. Assuming that charm contributions can still be 
neglected at our input scale, we adopt from global PDF fits a SU(3) fiavor asymmetric sea with 

5 

and we find so at the input scale the following SU(4) fiavor nonsinglet multiplets: 

(3) ^ ~ , hI'^ = H''^+^ + Hf+^ -2Hf^^ ^-Hr\ (24) 

J J J i/ J J J 



^(15) ^ Hf'-^ + Hf-'^ + Hf^^-3Hf^^^H]'\ 



which will be evolved autonomously, while as mentioned afore the singlet contribution Hj' ^ HJ^^ 
will mix with gluons. After evolution we finally return to individual quark fiavors. 

Hence, we have only to model the (conformal) moments Hj for the net sea quark and gluon 
contributions at the input scale. In both cases we take for the PDF Mellin moments a simple, 
however, realistic ansatz and we decorate it with t-dependency 

m,^')-N ''^l-"^'f^\'\!^'-\ ait), W = 0) = 1. (25) 

B[2 — a,p + l) l+j — a — a't 
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which for t = corresponds to a — x)^ PDF ansatz with momentum fraction average A^. 

Obviously, a and (3 determines the small and large x behavior, respectively. In accordance with 
phenomenological findings we fix the /3 parameters to be slightly larger as given by counting rules: 

(3'^'' = 8, and (3^ = 6. 

while a^'^^ ~ a*-^ > 1 is taken to be the effective intercept of the "pomeron". Note that the 
momentum sum rule implies the constraint 

In accordance with phenomenological findings from global PDF fits we set the averaged momentum 
fraction of u and d valence quarks to A^™^ = 0.4 and together with the intercepts for sea quarks 
and gluons, contained from a simple PDF fit [40j to HERA data, we fix the corresponding PDFs 
at the input scale Ql = 4GeV^: 

A^'^'^ = 0.152, a'^<'^ = 1.158 and A^^ = 0.448 , = 1.247. (26) 

The t-dependency of the GPD moment ( |25ll is contained in both the leading "Regge" (or "pomeron" ) 
pole l/(j + 1 — a — a't) and the residue 

/3(t)G{e^*,(l-t/m2)-'}, (27) 

chosen to be exponential with slope i? or as a dipole with cut-off mass m. The "pomeron" slope, 
observed in electroproduction processes, is smaller than the soft pomeron one a'p = 0.25/GeV^ 
and we set it here to 

^/sea _ ^/G _ o.l5/GeV2 . (28) 

The typical value of the slope parameter is for the processes of interest at the input scale Qq = 
4GeV^ measured to be -B = b/2 ~ 3 GeV^ and decreases with growing Q^. As noted above, its 
value for the gluon dominated exclusive J/ip electroproduction is considerable smaller B'^^'^ = 
~ 2.2GeVl In our fits we will prefer the dipole ansatz (1271) for the residue /3(t), where we 
might take the cut-off masses as in DVCS fits [40] 

M'^^ = ya5GeV, =\f^G(N. (29) 

3 Collinear factorization versus measurements 

To confront the collinear factorization approach with DVCS and DVMP data, we evaluate, as 
described in the previous Section, the differential cross sections ([1]) and (jl]) in terms of our GPD 
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moments fll5ll25p . The t-integrated cross sections are obtained from the differential ones by 

a{W, = t'^^'\t Q^) ^ , M = 0.938 GeV , (30) 



dt 



CBj- — 



where |tmin| ~ and |tcut| ^ 1- The bulk of DVMP data have been provided by the HI and ZEUS 
collaborations for the total cross sections 

cr = ax + ecTL, 

where the photon polarizability 

y= ^^^o , = 300 GeV [HERAII : 320 GeV] . (31) 

I — y + ^y^ s — 

To employ in our analysis this larger data set, we use the experimentally extracted R = ctl/ut 
raticl^, in form of a simple fit, shown in Fig. [H 

pexp.e2^ QV< ^ith /« = 2-2, P = 0.451, mv = 0.776 GeV 1 f P'\ (on) 

^ (1 + aQVO^ ia = 25.4,p = 0.180,mv = 1.02GeV / ^°M0J ^ 

where a possible weak W dependency and t dependency is neglected. Although the parameters a 
and p are strongly correlated, these fits indicate that higher twist-effects, parameterized as 

.2 \ p-i 



^exp 



are weaker suppressed than the canonical 1 / expectation [32] . One might imagine that this 
modification arises from resumed logarithmical corrections, which are expected from the break- 
down of factorization for the transverse polarized cross section [32]. Our predictions, e.g., for the 
t-integrated cross sections (jSOD are then obtained from (jl]) and (15^ . 



a{W, Q' 



e{W, Of) + 



'^^ ■ (33) 



In utilizing them, we will not take into account the errors from the i?-ratio fit fl32p and, moreover, 
as motivated in Sect. [T]we will ignore data points at lower photon virtualities, i.e., for < 4 GeV^. 



Since we replace here the flavor symmetric sea by the more realistic flavor scheme (123!) and we 
include recent DVCS data, we should first update previous DVCS fits [101111]. Thereby, the GPD 
parameters might be correlated, in particular, this is the case for the two skewness parameters 
and cut-off mass for both sea quarks and gluons. Note that evolution allows us to access partially 

^ In the ZEUS analyses the hypothesis of s-channel hehcity conservation was employed [71 HJ], while the HI 
collaboration used an improved approximation [5]. 
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Figure 1: i?-ratio from the HI (filled circles) [6] and ZEUS (filled triangles) [ZIIIT] collaborations for p*^ 
(left) and (j) (right) production, where statistical and systematical errors are added in quadrature. The 
solid curves show our fits (|3T]) . 



the gluon GPD, however, with large uncertainties. To reach convergency, we use as in previous 
DVCS fits |10] the PDF parameters (126|) . "pomeron" slopes ( 128|) . cut-off masses (1291) . and gluonic 
skewness parameter s'f from the KM 10b fit [H]. From a three parameter fit to the DVCS cross 
section (fT|5ni) measurements of the HI [TTl and ZEUS ITS] collaborations we find with 
X^/d.o.f. ~ 130/123 the new skewness parameters 

a'^"" = 1.158 , = -0.550 [-0.460] , = +0.130 [+0.094] , M''^ = VcTs GeV , (34) 

= 1.247, = -2.397 [-2.515] , = +0.892 [+0.892] , =V(l7GeV, 

which slightly differ from the KM 10b ones, given in square brackets. The DVCS fit is displayed in 
the three panels of Fig. [2] as dash-dotted curves, where the t-, W-, and Q^-dependencies are well 
described. Note that the mismatch between the dimensional counting prediction of a 1/ fall-off 
for fixed xb^, see cross section ([1]), and the measurements of roughly one power [161 [H] is resolved 
by the perturbative prediction of scaling violations. This prediction depends also on the chosen 
parameterization of the non-perturbative distributions at the input scale. Contrarily to PDFs, 
where the evolution at small x is essentially determined by the chosen value of the "pomeron" 
intercepts for gluons and the value of the input scale, the GPD evolution is also controlled by the 
values of skewness parameters S2 and S4 as well as to some extend by the different t-dependencies 
of flavor singlet quark and gluon GPDs. We also show predictions from the GK07 model [37] 
as dotted curves, where the GPDs were build from Radyushkin's double distribution ansatz [19] 
and adjusted to electroproduction data of light vector mesons. The world DVCS data set for 
small is well described by this x^/n.o.p ^ 226/126 prediction, where the model provides an 
almost perfect LO description and only fails to describe the W^-dependency (middle) of the low 
= 2.4 GeV^ ZEUS data (filled circles). 



9 




Figure 2: Differential DVCS cross section vs. —t (left) as well as t-integrated ones vs. W (middle) and 
vs. (right) are taken from [16] (filled rhombus), [T7| (empty squares), [18] (filled circle), and [19] 
(filled squares), where statistical and systematical errors are added in quadrature and normalization 
uncertainties were ignored. Measurements are confronted with a DVCS fit (j34p (dash-dotted), the GK07 
model prediction (dotted) [37|, and simultaneous DVCS/DVMP fits ([35]) (dashed) and ([36]) (sohd). 

A large set of DVMP data at small xbj were obtained in the p° channel from the HI [H |3l [5l [6] 
and ZEUS [2l HI |7] collaborations (0 measurements are listed in [9l El [10] and [H IH]). The most 
recent publications from ZEUS [7j and HI [0] refer to p° data from the 1996/97 and 1998-2000 
HERA runs. The t-integrated cross section measurements for fixed W = 75(90) GeV versus Q^, 
where the upper t-cuts are = 1 GeV^ and |)f:^\| = 0.5 GeV^, are consistent with each 

other, see filled up-triangles [6] and empty down-triangles [7] on the right upper panel in Fig. [31 
The measurement of the longitudinal cross section (empty up-triangles) has been achieved by 
the HI collaboration from the knowledge of the spin density matrix elements, shown on the 
right p° and (p panels of Fig. [31 Thereby, the experimental errors slightly increase due to the 
uncertainties of the i?-ratio, see Fig. [H One also realizes from these panels that our DVCS 
predictions (dash-dotted curves) overshoot the DVMP cross sections and falls off too steeply 
with growing Q^. Experimentally, the power-like fall-off of the cross sections is determined as 
~ 1 / [6] while dimensional counting predicts a 1/ fall-off, see the perturbative prediction 
(Hl[5l) . However, it might be too naively to conclude that these discrepancies already rule out the 
collinear factorization approach rather they might be attributed to our relative hard gluon GPD 

Compared to the HI measurements [6], both the t- and TV-dependency of the ZEUS measure- 
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Figure 3: Differential cross section versus —t (left) as well as t-integrated ones versus W (middle) 
and (right) for DVMP of (up) and (p (down). Longitudinal (empty up-triangles) and total cross 
sections (filled up-triangles) for and cj) meson production form HI [B] and ZEUS (empty down-triangles) 
[71 [11]. Statistical and systematical errors are added in quadrature and normalization uncertainties 
were ignored. Measurements are confronted with a DVCS fit (j34p prediction (dash-dotted) and two 
simultaneous DVCS/DVMP fits shown as dashed (p5|) and solid (j36|) curves. 

ments [7] are slightly flatter, e.g., 

r b^\Q' = 11.5 GeV^) \ ^ r 6.72 ± 0.53 \ 
j5ZEUS(g2 = nGeV^)J 15.7 ±0.5 l":^ // ' 

r 5Hi(g2 ^ g g Ggv^) ^ J 0.57 ± 0.10 1 

1 52EUS(g2 = g GeV') J ~ I 0.4 ± 0.052 ^^0!! J ' 

see also filled up-triangles and empty down-triangles on the upper left and middle panels in Fig. [3J 
A slightly flatter ly-dependency of the ZEUS data is also established in the (p channel [11], see 
lower middle panel in Fig. |3l These differences might be attributed to systematic uncertainties 
of the background subtractions, in particular of proton dissociative contributions that were ex- 
perimentally studied by the HI collaboration [5]. Although HI and ZEUS data are compatible to 
each other, their separate uses imply some freedom in the partonic interpretation. 
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The DVMP data from ZEUS [TJ [TT] are describable to LO accuracy with a very soft gluon 
GPD, where its leading "Regge" intercept at the input scale Qq = 4 GeV^ is even smaller than one, 
i.e., < 1. This ensures that the increase of the effective "Regge" intercept with growing Q^, 
which is driven by the j = pole of the anomalous dimension in the gluon channel, is sufficiently 
slow. Furthermore, it turns out that such a very soft gluon GPD is also compatible with the 
DVMP measurements from HI and the DVCS data set. This is illustrated by the dashed curves 
in Figs. Eland E], which arise from a simultaneous DVCS/DVMP fit with x^/d.o.f. ^ 618/297 that 
pins down the eight parameters 

^sea _ ^_^g^ ^ ^sea _ Q ggg ^ ^sea _ _Q 21Q , M^"^ = 7(^554 GcV , (35) 

=0.513, = 1.950, =-0.469, = VO.462 GeV . 

The new sea quark intercept and cut-off mass are entirely consistent with the DVCS fit (!34l) . 
The skewness parameters have now a reversed sign, providing us an alternative solution to 
the optimization problem. However, the very low gluonic intercept might be inconsistent with 
PDF findings from deep inelastic scattering (DIS) measurements. If we fix this intercept, e.g., 
= 4GeV^) = 1.1, the rather fiat VT-dependency of the ZEUS data implies a mismatch in 
the normalization. Hence, in such combined HI and ZEUS DVMP fits we can only get disfavored 
X^/d.o.f. ~ 5 values, see also the pomeron fits in [50j. 

Let us now rely on the DVMP measurements of the HI collaboration [6J and DVCS data 
[lEl [ni [TBI [19]. If we assume a softer gluon PDF with a^{Q'^ = AGeV^) = 1.1, as it, e.g., also 
appears in the hand bag model [37], we can reach a good simultaneous DVMP/DVCS description. 
Thereby, we adopt the quark PDF from the DVCS fit flMj) and ask for the remaining six skewness 
and dipole cut-off parameters. From a x^/d.o.f. ~ 342/224 fit, shown as solid curves in Figs. |2] 
and [3], we find 

a'"'' = 1.158 , = +0.802 , = -0.278 , M'"'' = VO.548 GeV , (36) 

=1.100, =-4.288, =+1.616, = VO.351 GeV. 

Since the three gluonic parameters and M*^ are strongly correlated, the t-dependency of 

the gluon GPD can even in this simultaneous DVCS/DVMP fit not be pinned down. 

Let us also mention that the u channel might be reasonable described by our two simultaneous 
DVCS/DVMP fits. For both of them we find that the longitudinal cross section ratios u)/p^ ~ 0.1 
at = 70 GeV are compatible with the measured ones [T2|, e.g., 

_ (W = 70 GeV, = 7GeV^) = 0.089 + 0.014 + 0.019. 

In the left panels of Fig. H] we compare our models with a standard PDF parameterization of 
Alekhin [51] . Our total sea quark PDF from the DVCS fit (dash-dotted curves) is the same as in 
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Figure 4: PDFs (left), GPDs on the cross-over line (middle), and the skewness ratios (j37p for the flavor 
singlet sea quark (up) and gluon (down) models, employed in Figs. [2] and [3l Phenomenological PDFs at 
LO (grayed area) and NLO (light grayed area) are taken from Ref. |51j . 

the simultaneous DVCS/DVMP fit ( l36l) and it is compatible with Alekhin's LO parameterization 
(grayed area). To some extend this is also the case for the other sea quark models (dashed and 
dotted curves). Note also that the NLO radiative corrections in the sea quark sector induce only 
mild reparameterization effects, compare grayed and light grayed error bands. The sea quark 
GPDs on the cross-over line, shown in the upper middle panel, are consistent with each other and 
so we might conclude that our sea quark models are constrained by DVCS data. Although the 
skewness parameters in the DVCS fl34|) and simultaneous DVCS/DVMP f p6|) fits have different sign 
and magnitude, the resulting sea quark GPDs, shown as dash-dotted and solid curves, are hardly 
to distinguish at the input scale. Also the gluon PDF in our DVCS fit (dash-dotted curves) is 
compatible with standard LO parameterizations. The gluon PDF of the GK07 model [37] (dotted 
curves) and in our simultaneous DVCS/DVMP fit (136|1 are rather similar and underestimate the 
phenomenological LO findings. It is worth to mention that radiative corrections will drastically 
reduce the gluon PDF and, hence, these both aforementioned models are compatible with NLO 
PDF parameterizations. The gluon GPD on the cross-over line is for the GK07 model roughly 
given by the PDF, while in our more flexible models the GPDs on the cross-over line are much 
smaller than the gluon PDF. Our simultaneous fit to DVCS and Hl/ZEUS DVMP data results in 
a very soft gluon PDF (dashed curves), which is inconsistent with phenomenological PDF findings. 
In the right upper and lower panels of Fig. H] we display the quark and gluon skewness ratios 

and ,37) 
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for fixed x = 10'^ versus For our DVCS §S) and simultaneous DVCS/DVMP ([Ml) fit we 
find as previously the value r^'^'^ 1 [SSI- The sea quarks from these both models, which mostly 
coincide at the input scale, evolve only slightly. The stability of the sea quark ratio under evolution 
requires that the corresponding gluonic r-ratios are smaller than one. In the coUinear factorization 
framework to LO accuracy this GPD feature is needed for a successful DVCS fit [52j. The sea 
quark ratio of the GK07 model (dotted curves) is essentially larger, r^°^ ^1.3 and rather stable 
under evolution, too. Note that the successful DVCS description of the GK07 model presumably 
originates from the interchange of evolution and skewing procedure [37], see numerical examples 
in [53] and comments in Sect. 3.1 of [40] . 

4 Summary and conclusions 

Based on the coUinear factorization approach at LO accuracy and flexible GPD models, we demon- 
strated that simultaneous GPD fits with x^/d.o.f. ~ 2 or better describe the kinematical variable 
dependency of DVMP (light vector mesons) and DVCS measurements in the small xbj region 
already for a photon virtuality of > 4GeV^. In our studies we were left with some theoretical 
and experimental uncertainties. For instance, we did not extensively explore the different partonic 
degrees of freedom that enter in these processes, in particular, we fixed the flavor content of the 
quark sea and we relayed on the asymptotic form of DAs, which do not change under evolution. 
Moreover, we used naive renormalization and factorization scale setting prescriptions. To utilize 
the larger data set for the unpolarized DVMP cross sections, we used the experimental i?-ratio, 
where the hypothesis of s-channel helicity conservation was employed and it was assumed that 
this ratio only depends on Q^. Furthermore, we simplified our analysis by ignoring errors in the 
i?-ratio and normalization uncertainties in the data sets. 

In our studies we did not encounter difficulties in the unifying description of t-dependency, how- 
ever, the inclusion of DVMP data from the ZEUS collaboration challenges the expected "pomeron" 
like behavior of the gluon GPD. The successful description of these data requires a very soft gluon 
GPD at the input scale of 4 GeV^, which might be ruled out by standard PDF parameterizations. 
Employing only the DVMP measurements from the HI collaboration together with the DVCS 
data allow us to describe these exclusive channels, however, with a rather soft gluon GPD that on 
the cross-over line is smaller than the gluon PDF. Such a skewness effect at LO accuracy has been 
already observed in DVCS fits with more realistic gluon PDFs [iQl HI]. Hence, we expect that 
a reasonable global description of the full DIS, DVMP and DVCS data set at small xbj cannot 
be reached or is disfavored at LO level. On the other hand if we would restrict ourselves to the 
few released HI data points for the longitudinally t-integrated cross sections, this task might be 
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succeeded. Nevertheless, from our analyses we might conclude that t-channel quark exchanges 
are important in the DVMP processes. While the t-dependency of the sea quark GPD is rather 
well constrained from DVCS data, we certainly realized that the gluon GPD suffers from large 
uncertainties. 

Our partonic interpretation of DVMP differs from those of color dipole models in which by 
definition the utilized "microscope" is tuned to the gluonic component of the nucleon. Note that 
the separation of flavor singlet quark and gluon degrees of freedom is implicitly done by adopting 
a factorization scheme. In the handbag model approach higher twist contributions, in terms 
of transverse degrees of freedom, allows to adopt the popular Radyushkin's double distribution 
ansatz for the description of experimental DVMP data. The resulting GK07 GPD model is 
qualitatively different from our ones. Nevertheless, it also provides a good LO description of 
DVCS measurements, which originates from the specific modeling of GPD evolution. 

Certainly, in our partonic description we are left with some discrepancies between DIS, DVCS, 
and DVMP findings, which in spite of experimental, theoretical, and model uncertainties cannot be 
taken as a convincing counter argument against the coUinear factorization approach. In our opinion 
it is worth to study such a global fitting procedure in the NLO approximation of this approach. 
Thereby, one should also include electroproduction data of J/ip measurements which strongly 
constrain the gluon GPD. The software tool that is needed for this task is under development 
and, certainly, the fitting procedure should be improved by taking into account the uncertainties 
of both the ai / ct separation and the cross section normalization. 
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